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A candidate for a consistent steady state thermodynamics is constructed for a radiation field 
in vacuum sandwiched by two black bodies of different temperatures. Because of the collisionless 
nature of photons, a steady state of a radiation field is completely determined by the temperatures 
of the two black bodies. Then the zeroth, first, second and third laws can be extended to steady 
states, where the idea of local steady states plays an important role for the system whose geometrical 
shape is anisotropic and inhomogeneous. The thermodynamic formalism presented in this paper 
does not include an energy flux as a state variable. This is consistent with the notable conclusion 
by C. Essex, Adv. Thermodyn. 3 (1990) 435; Planet. Space. Sci. 32 (1984) 1035 that, contrary to 
the success in the irreversible thermodynamics for dissipative systems, a nonequilibrium radiation 
field does not obey the bilinear formalism of the entropy production rate using an energy flux and 
its conjugate force. Though the formalism given in this paper may be unique to a radiation field, a 
nonequilibrium order parameter of steady states of a radiation field is explicitly defined. This order 
parameter denotes that the geometrical shape of the system determines how a steady state is far 
from an equilibrium. The higher the geometrical symmetry, the more distant the steady state. 
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I. INTRODUCTION 



The traditional treatment of the radiative energy transfer [l| has been applied to a mixture of a radiation field 
with a matter like a dense gas or other medium. In such a traditional treatment, the successive absorption and 
emission of photons by components of the medium matter makes it possible to consider that the photons are as if in 
local equilibrium states whose temperatures equal those of local equilibrium states of the medium matter. However, 
because photons are collisionless, it is impossible to apply the idea of local equilibrium to a radiation field in vacuum. 

By the way, when one deals with nonequilibrium states for a matter (without including a radiation field) , it is usual 
to place an interest on the heat flux. Indeed for example, the classical and extended irreversible thermodynamics 
Q make use of the bilinear formalism of the entropy production rate using the heat flux and its conjugate 
thermodynamic force to obtain many developments in describing nonequilibrium phenomena. Then consider the case 
that a matter is put into a space where nothing but a nonequilibrium radiation field exists. This system possesses 
an interaction of the matter with the outside radiation field by the absorption and emission of photons at the surface 
of the matter. Note that, concerning a thermodynamic treatment of this system, it has already been revealed in 
references [BJ that, because of the collisionless nature of photons in vacuum, the entropy production rate for the 
whole system composed of the matter and the radiation field can not be expressed in the bilinear form. And it is 
also stated in Q [j| that, even if the bilinear formalism was applied to this system, disappearing radiative energy 
flux (or thermodynamic force of the radiation field) would not always denote disappearing thermodynamic force (or 
radiative energy flux). This means that the energy flux does not work as a consistent state variable for a system 
including a nonequilibrium radiation field in vacuum. Then two conclusions follow, that the traditional treatment 
of the radiative transfer mentioned above is applicable only to a mixture of a radiation field with a matter which is 
dense enough to ignore the vacuum region among components of the matter, and that a distinct research from that of 
ordinary dissipative matter systems is required to understand a nonequilibrium radiation field in vacuum. However, 
a consistent thermodynamic formulation for a system including a nonequilibrium radiation field in vacuum has not 
been accomplished, and a consistent nonequilibrium order parameter for a radiation field has not been obtained so 
far. 

This paper shows a candidate for a consistent steady state thermodynamics for a radiation field, and obtains the 
consistent nonequilibrium order parameters of the steady states for a radiation field. We consider a radiation 
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field in vacuum sandwiched by two black bodies of different temperatures. Keeping each temperature constant, the 
radiation field is in a nonequilibrium steady state with a stationary energy flow given by the Stefan-Boltzmann law. 
The settings treated in this paper are illustrated in figure Q Here it should be emphasised that, as pointed out at §63 
in reference 0, the collisionless nature of photons denotes that a nonequilibrium radiation field can never relax to 
any equilibrium state unless there exists any interaction with other matter. Therefore the existence of the two black 
bodies in our settings is essential in formulating the steady state thermodynamics for a radiation field. Especially 
in extending the second law to steady states, it is necessary to consider a relaxation process of the whole system 
composed of the radiation field and the two black bodies. Given a well defined steady state entropy of the radiation 
field, the total entropy of the whole system have to be monotone increasing during the relaxation process of the whole 
system. 

Before proceeding to the construction of the steady state thermodynamics for a radiation field, it is helpful for 
readers the steady states of ordinary dissipative systems and the fluctuation of our radiation system are mentioned. 
As mentioned above in second paragraph, the thermodynamic formulation presented in this paper is unique to a 
radiation field. Therefore, when one attempts to construct a steady state thermodynamics for dissipative systems 
(without including a radiation field), a consideration different from ours should be made. Indeed, the reference Q 
suggests a consistent macroscopic formulation of steady state thermodynamics for some dissipative systems like heat 
conduction, shear flow, electrical conduction and so on. In their formalism, the heat flux plays an essential role 
as a consistent state variable contrary to our radiation system, although it is also reported |8( that their steady 
state thermodynamics is not omnipotent and is applicable to some restricted class of dissipative systems. Further in 
reference , the velocity distribution function of the steady state Boltzmann equation for a hard core molecule gas 
has already been explicitly derived up to second order in the density and temperature gradient. Their distribution 
function does also depend on the heat flux. Here we should note that, although it may appear strange that our 
formalism for a radiation field does not include an energy flux as a state variable, it seems that the heat flux plays 
the role of a consistent state variable for dissipative systems. The absence of an energy flux from consistent state 
variables should be understood as a unique property to a radiation field. 

We next turn to the fluctuation of our radiation system. Because photons are collisionless, the distribution functions 
of photons in our settings are constructed by simply modifying an equilibrium Planckian distribution as shown in 
following sections. Consequently, the steady state of the radiation field is completely determined by the temperatures 
of the two black bodies. In addition, the collisionless nature shows that no fluctuation arises from the radiation 
field itself. The only possibility of a fluctuation arises from the absorption and emission processes of photons at the 
surfaces of the two black bodies. If the energy flux in the radiation field becomes strong enough and if the time 
scales of thermalisation of two black bodies are long enough, then the black bodies may be affected as to go into 
nonequilibrium states. In such a case, nonequilibrium temperatures of the bodies are unknown in general, and the 
distribution function of photons should be determined according to nonequilibrium statistical prop erties of the bodies. 
Concerning such a totally nonequilibrium case, there have already been some reports [Tl| |12| trying to determine 
the distribution function for a radiation field emitted by a nonequilibrium matter. They consider a radiative energy 
transfer inside a dense matter, and apply the information theory to the whole system composed of a radiation field 
and a medium matter. pol | However, after those reports were published, it has been revealed in reference |l3j that, 
at least for a matter whose components are colliding and interacting with each other, the distribution function for a 
steady state of the matter derived with the information theory does not qual itatively agree with that derived with 
a steady state Boltzmann equation. Further it has also been concluded in [13j that the nonequilibrium temperature 
determined with the information theory has no physical meaning. Therefore, because their distribution function of 
the radiation field depends on the nonequilibrium temperature of the medium matter derived with the information 
theory, the reliability of the distribution function of the radiation field may not be given in those reports 0] 01 ■ 
Hence there is no confirmed form of the distribution function of a radiation field emitted by a nonequilibrium matter. 
In this paper, as a first step to construct a consistent macroscopic theory for a radiation field, we simply assume that 
each black body is always in an equilibrium state and no fluctuation arises. 

The contents of this paper are organised as follows. In section \n\ our basic requirements are explained. Then 
the steady state thermodynamics is established for the cylinder case (left one in figure ^) in section IIIII where a 
nonequilibrium order parameter of steady states is defined. Section IIVI is devoted to the extension of the results 
obtained in section ITTT1 to the cavity case (right one in figure^), where the idea of local steady states plays an essential 
role. Summary and discussions are in section where special properties of a radiation field, the distance of steady 
states from equilibrium states, future improvements and possible applications are discussed. 

Throughout this paper, c is the speed of light, h = h/2n is the Planck constant, a — ir 2 /60h 3 c 2 is the Stefan- 
Boltzmann constant and the Boltzmann constant unity ks — 1 is set. 
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Cylinder case Cavity case 



FIG. 1: A steady state for a radiation field. J is a stationary energy flow. 



II. BASIC REQUIREMENTS AND CONDITIONS 

We begin with the following basic requirements for constructing the steady state thermodynamics (SST) for a 
radiation field: 

(Rl) The state variables of equilibrium thermodynamics are extended to steady states. And at equilibrium limits, 
those extended variables take on the same value as the equilibrium ones. 

(R2) There exists a state variable which describes a nonequilibrium order of steady states. And that variable 
disappears at equilibrium limits. 

These requirements denote the zeroth law of the SST, the existence of steady states. A steady state of a radiation 
field is realised in the settings shown in figure In the cylinder case, the left and right lids are black bodies and in 
equilibrium individually with constant temperatures Tl and Tr respectively, where the side wall of the cylinder is a 
perfect mirror of zero temperature and no heat is exchanged between the side wall and each lid. Here the side wall is 
considered just for preparing a radiation field of finite volume. Therefore we ignore any thermodynamic effect of the 
side wall hereafter. In the cavity case, the inner and outer black bodies are in equilibrium individually with constant 
temperatures Ti n and T ou t respectively. Consequently due to the Stefan-Boltzmann law, a stationary energy flow J 
exists in both cases. The equilibrium limit is taken with Tl = Tr and T jn = T out . 

The requirement of the zeroth law is not enough to search for consistent definitions of state variables. Therefore 
we require that the following properties of equilibrium thermodynamics are retained in the SST for a radiation field: 

(R3) Any state variable is either intensive or extensive. 

(R4) The free energy is convex (or concave) with extensive (or intensive) variables. 
(R5) Thermodynamic functions are related by the Legendre transformation. 
(R6) The Gibbs-Duhem relation holds to ensure the first law. 
(R7) Extensive variables are additive. 

To make the meaning of requirement (R7) more clear, consider the cylinder case as an example. The internal energy 
of the left lid, El, and that of the right lid, En, are determined by equilibrium thermodynamics, since the lids are 
in equilibrium individually. Then with the SST internal energy, E ra d, of a radiation field defined in section ITTT1 the 
requirement (R7) denotes that the total energy is given by E to t = El + Er + E rac i- Similarly, the total entropy is 
Stot = Sl + Sr + S ra d, where Sl and Sr are equilibrium entropies of left and right lids, and S ra d is the SST entropy 
defined in section lTTTl We require that the same additivity holds for the cavity case. After introducing the idea of local 
steady states for the cavity case, the additivity of extensive variables of the radiation field in different local steady 
states is also considered in section llVl 

Though the additivity for extensive variables of the radiation field is an appropriate requirement due to the colli- 
sionless nature of photons, it should be clarified that the additivity among extensive variables of the two black bodies 
and the radiation field like Stot — Sl + Sr + S ra d is a less appropriate requirement and is nothing but a simple 
assumption. However, as mentioned at the end of fifth paragraph in sectionQJ we assume the following condition: 

(CI) Each black body is always in an equilibrium state, and no fluctuation arises on the radiation field. 

Under this restriction, the state variables of the two black bodies are always given by the ordinary thermodynamics 
(whose extensive variables are additive), and it is appropriate to consider that, once a photon is emitted at a black 
body, no interaction of the photon with black bodies exists until it is absorbed at a black body. This implies that 
the interaction between the radiation field and a black body is of a very short range, and that it seems not so bad to 
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FIG. 2: Relaxation process. Total entropy should be non-decreasing, dStot > 0. 

assume the additivity like Stot — Sl + Sr + S ra d- When one proceeds to the study of the fluctuations in our settings 
which occur in the case when the bodies are in noncquilibrium states, the possibility of existence and influence 
of non-additive state variables of the bodies should be reconsidered. And when the bodies are so massive that the 
gravitational effect which is a long range interaction between the bodies can not be ignored, the possibility of existence 
and influence of non-additive state variables of the bodies should be reconsidered. Further if the gravity becomes so 
strong that the general relativistic effects can not be ignored, a non-additivity among state variables of the bodies 
and of the radiation field may arise. And concerning the condition (CI), a practical issue on future improvement is 
mentioned in section Ivl 

Next in order to extend the second law to steady states, recall the third paragraph in section HJ The radiation 
field in a steady state can never relax to an equilibrium state without introducing interactions with other matters 
. Therefore it is necessary for the second law to consider the relaxation process of the whole system composed of 
the radiation field and the two black bodies. Here as an example, let us consider again the cylinder case. Cover the 
cylinder with a heat insulator to be isolated. Then the whole system composed of the two lids and the radiation field 
relaxes to a total equilibrium state in which the lids and the radiation field have the same equilibrium temperature 
(figure [21 . If each lid passes a sequence of equilibrium states during this relaxation process, the radiation field 
passes a sequence of steady states. The same relaxation process of the whole system in the cavity case is also easily 
considered. As mentioned at the third paragraph in section [H given a well defined SST entropy S ra d, the total 
entropy Stot (= Sl + <Sr + S ra d) have to be monotone increasing during such a relaxation process of the whole system. 
Therefore we have to place the other requirement: 

(R8) Total entropy is non-decreasing, dStot > 0, along the relaxation process of the whole system. And the equality, 
dStot — 0, holds for a total equilibrium state. 

Though it may appear strange that the non-decreasing nature of the entropy of a nonequilibrium radiation field is not 
given with a relaxation of an isolated radiation system, such a property should be understood as a unique property 
to a radiation field. Further, in order to realise the relaxation process of the whole system in a simple way, we require 
the following condition: 

(C2) The time scale for changing the equilibrium states of the two black bodies is sufficiently longer than the flight 
time of a photon in the space between the two black bodies. 

If this condition is not satisfied, the temperatures of the black bodies may change while a photon travels in the space 
between the two bodies. That is, the effects of retarded time on photons should be taken into account. However the 
condition (C2) allows us to ignore the retarded time effects. 

Further concerning the relaxation process of the whole system, we should pay attention to the special relativistic 
effect as follows. If the volumes of the two black bodies change during the relaxation process, their surfaces may 
move. This denotes that the photons are emitted from a moving source and that the special relativistic Doppler effect 
due to the Lorentz transformation should be considered. However, because the coefficient of thermal expansion of 
an ordinary solid is small in general, the volumes of the two black bodies do not change so much and the speeds of 
their surfaces are so small in comparison with the speed of light that the Doppler effect is negligible. Therefore, for 
simplicity, we assume that the two black bodies satisfy the following condition: 

(C3) The two black bodies are rigid so that their volumes do not change under the influence of the energy flux in 
the radiation field. 

In summary, we aim to construct the steady state thermodynamics for a radiation field in the settings satisfying 
the conditions (CI), (C2) and (C3) under the requirements from (Rl) to (R8). 
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III. CYLINDER CASE 



A. Energy flux and scale change 



The cylinder case is treated in this section. A steady state of a radiation field is accompanied by a stationary energy 
flow determined by the Stefan-Boltzmann law, J = a ( — T# ) A, where A is the cross section of the cylinder. 
Hereafter we set Tl > Tr without loss of generality. 

As seen below, the steady state of a radiation field is completely determined by temperatures Tl and Tr without 
respect to the cross section A and the length of cylinder L. This motivates us to restrict the scale change to keep 
the cylindrical symmetry of the system. That is, the scale change is operated by changing the volume V (— AL) as 
V — > XV, where A > and the scalings of A and L are not uniquely determined. The extensivity and intensivity are 
defined with this scale change. Then, the volume of the radiation field, V, is obviously an extensive variable. However 
the energy flow J is neither intensive nor extensive, since J oc A but J 9^ L. Therefore J is not a state variable due to 
the requirement (R3). On the other hand, when one introduces an energy flux (an energy flow per unit cross-section), 
j = a [Tl — T£ ), this can be a candidate for a state variable. However it will be shown in subsection IIII II that j 
does not work as a consistent state variable. 



B. Distribution function 



Photons at the moment of emission at a lid have a Planckian distribution of temperature Tl or Tr, since the lids 
of cylinder are in equilibrium due to the condition (CI). Then, because photons are collisionless, the photons emitted 
at one lid do not alter their Planckian distribution until absorbed at another lid. Hence we can define a distribution 
function, d(p), of a steady state for a radiation field as 

55 cxp[^/r (p )] - 1 ' (1) 

where p is a momentum of a photon related to its frequency to — pc/h and T(p) is given by 

m = 



T L for p = pl 
T R for p = pn. 



where pl is the momentum of a photon emitted at the left lid and pu is that emitted at the right lid (see figure |3j). 

This d(p) implies that, under the condition (CI), the steady state of a radiation field can be considered as a 
superposition of two equilibrium states of temperatures Tl and Tr. Then one may think that there is no meaning in 
study on the radiation field under the condition (CI). However, note that the superposed radiation field as a whole is 
in a steady state with the energy flux j, and that it has already been revealed that j can not work as a state variable 
Therefore, the study of a radiation field under the condition (CI) obviously has two meanings to reveal an 
appropriate state variable of a nonequilibrium order, and to show an example of a consistent SST, although the new 
steady state variable and the way of construction of the SST presented in this paper are unique to a radiation field. 



C. Internal energy 



Wc define an SST internal energy of a radiation field, E rac [, as the energy carried by all photons, 

~ 3 h 



2 J ^-dx 3 hud(p) 



47T 3 C 3 



47T 3 C 3 



V 



dcjdQpLU 3 d{pL) 



duidQpdx 3 to 3 d(p) 
dujdQpLu 3 d(pn) 



where the factor 2 in the first line is due to the helicity of a photon, dVl p is the solid-angle element in p-space, and 
a relation J dp 3 /h 3 = (l/87r 3 c 3 ) j duidflpUi 2 is used. Since the integrands in the last line do not depend on fl p , the 

angular integrals give half of a total solid-angle, J~ dil p — dVt p = 2ir J* d9 sin 6 = 2ir, where 9 is an angle shown 
in figure [3] Then the SST internal energy E ra d becomes 
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V 



(E eq (T L ,V)+E eq (T R ,V)) , 



(2) 
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FIG. 3: Photons of momentum pl and pr. 



where J °° dx x 3 / (e x — 1) = 7r 4 /15 is used, and E eq (T, V) — (4<r/c) T 4 is the internal energy of a radiation field in 
an equilibrium state of temperature T and volume V. This E rad is obviously extensive and satisfies the requirement 
(Rl), E rad = E eq (T eq , V) for T L =T R = T eq . 



D. Entropy 

To search for a definition of an SST entropy, we try referring ^ij 0] and §55 in |(|. In those references, a general 
form of an entropy of a boson gas in a nonequilibrium state has been derived by making use of the Boltzmann principle. 
That is, by counting the number of states of bosons and using the Stirling's formula, the entropy of a nonequilibrium 
boson gas, Sboson, has been obtained 

/dp 3 
-foT dx3 9p,x [ ( 1 + Np t g ) In ( 1 + Npg ) - Nf t $ In Np t3 ] , 

where gp t $ and Npg are respectively the number of states and the average number of bosons at a point (p, x) in the 
phase space of the boson gas. It has also been shown in reference that the maximisation of Sboson (S Sboson = 0) 
for an isolated system give s the equilibrium Bose distribution function. 

In references jy] 0] |l5j| . concrete forms of gp$ and Npg are not specified, since an arbitrary system is considered. 
However for our system, they are determined; gp t s = 2 due to the helicity of a photon and Np t $ — d(p). Hence, we 
define the SST entropy of a radiation field, S rad , as 

S rad = 2 J ^ dx 3 [ ( 1 + d{p) ) In ( 1 + dip) ) - d{p) \nd{p)] 

= ^(T 3 +T 3 )V=l(S eq (T L: V) + S eq (T R ,V)), (3) 

where S eq (T,V) — (16(t/3c) T 3 V is the equilibrium entropy, and the following relations are used, 
f^dx [x 2 /(e x - 1)] ln(e x -l) = ll7r 4 /180 and J °° dx [x 2 /{\ - e~ x ) ] ln(l - er x ) = -tt 4 /36. This S rad is obvi- 
ously extensive and satisfies the requirement (Rl), S ra d — S eq (T eq , V) for Tl = T R = T eq . 

As discussed in section [H] a well defined entropy has to satisfy the requirement (R8). Indeed as shown in next 
subsection, S ra d of equation © satisfies (R8). Therefore S rad of © is a candidate for a well defined SST entror 
of a radiation field. However, there are many other forms of S ra d which can satisfy (R8). Further the reference 
conjectures well the inappropriateness of the Shannon-type form like Sboson as a nonequilibrium entropy at least for 
dissipative systems, and hence there is a possibility that S ra d of J2J is also not appropriate for the SST entropy of a 
radiation field (collisionless system). What we can say from above is that, while S ra d of (JHJ satisfies the requirement 
(R8), it is not clear at present whether or not S rad of 10 is unique and completely appropriate to the SST entropy 
of a radiation field. Hence, to be exact, we observe whether we can construct a candidate for a consistent SST for a 
radiation field with S rad of J3J . 



E. 2nd law with a relaxation process 

This subsection shows that S rad of equation satisfies the requirement (R8). We set the cylinder being covered 
with a heat insulator and consider that each lid passes a sequence of equilibrium states and the radiation field passes 
a sequence of steady states during the relaxation process of the whole system. Due to the condition (C3), the volumes 
of the lids does not change during the relaxation. Then the following equilibrium thermodynamic relations hold for 
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FIG. 4: The case of dV 7^ during the relaxation process. The region denoted as EX is of true vacuum of zero temperature, 
where even a radiation field does not exist. That is, the smoothly movable wall never give work on the region EX. 



the lids, 

dE L = C L dT L = T L dS L , dE R = C R dT R = T R dS R , 

where El, Cl and Sl are the equilibrium internal energy, heat capacity and entropy of the left lid, and E R , C R and 
S R are those of the right lid. Since the whole system is isolated by the heat insulator, the total energy conservation 
holds, 

Etot = E L + E R + E rac t = const. 
=> = C L dT L + C R dT R + ^ ( TldT L + T*dT R ) V + ^ ( T L 4 + T* ) dV . 

Here note that, even if the volumes of the lids do not change due to the condition (C3), there is a possibility of the 
change of the radiation volume V during the relaxation process. Two examples of such a case are shown in figure 0] 
where the smoothly movable wall is pushed by the radiation pressure. That is, dV > during the relaxation. And 
it is appropriate to assume that the speed of the wall during the relaxation is small enough to ignore the special 
relativistic Doppler effect on photons and the kinetic energy of the lids. However it should be emphasised that, 
when the volume V continues to expand, the whole system never reaches any equilibrium state, since the equilibrium 
state does not evolve in time. Furthermore if the exterior region denoted as EX in figure 0] is sufficiently large, the 
radiation volume V becomes so large that the condition (C2) will be violated after a sufficiently long time passes. 
Hence it is appropriate to take the equilibrium limit with Tj, = T R and dV = for the relaxation of the whole system. 
Consequently, the differential of the total entropy Stot — Sl + S R + S ra d becomes 

dStot = ^dTL + ^dT R + —(TL 2 dTL+ndT R )V+^(T L "+n)dV 

1 L 1 R C o C 

The inequality dS to t > holds obviously for Tl > T R (dT R > 0), and with the help of the energy conservation it is 
also easily derived that the inequality dStot > holds for the case T R > Tl (dTL > 0). Further the equality, dStot = 0, 
holds for the equilibrium limit Tl = T R and dV = 0. As mentioned at the requirement (R8) in section^] this denotes 
that S ra d of equation is a candidate for a well defined SST entropy. 
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7 T r , T~l 



Cr 



dTn 



2<T 



T R d V 



F. Pressure 



We define an SST pressure by the tensor, 

P^ = 2j^-^ Cfl c u d(f), (4) 

where the factor 2 is due to the helicity of a photon, c M is a spatial component of a photon velocity satisfying 
E^c^ = c 2 , and fi and v denote the Cartesian components. This tensor P corresponds to the pressure tensor of an 
ordinary dissipative system, since p/c corresponds to the mass of a molecule. P is obviously symmetric, P M „ = P Vfl . 
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At the equilibrium limit Tl = Tr = T eg , the SST pressure tensor becomes 

P-itrP^ , trP = 2/ * _ , . g T,t = 3P.,(T W ) , 

where pc = fiw is used, trP is the trace of P and P eq (T) = (4cr/3c)T 4 is the equilibrium pressure. This denotes the 
SST pressure tensor P satisfies the requirement (Rl). 



G. Free energy 



We try to define an SST free energy, F rac i, as a scalar variable satisfying the following equilibrium limit, 

dF ra d 



dV 



Peq (T et 



Therefore, since V is also a scalar variable, we are to require 

dFrad 



dV 



= Prad , (5) 



where P ra d is the scalar quantity which is constructed using only the SST pressure tensor P and reduces to P eq at 
equilibrium limits. This scalar P ra d can be given by the trace of P as 

Prad = I trP = § / fc, d{p) = \ = I? (Ti + n) . (6) 



3 3 7 /1 3 K1 ' 3 V 3 c 

It is obvious that P ra d = Peq(T eq ) for Tl = Tr = T eq . As will be seen later, this definition of the scalar P ra d ensures 
the consistency of our SST formulation. Hence we define F ra d as 

Frad = -\ trP V = ~ ( Tt + T* ) V = \ ( F eq (T L , V) + F eq (T R , V) ) , (7) 

where F eq (T,V) — — (4a /3c) T 4 V is the equilibrium free energy. This F ra d is obviously extensive and satisfies the 
requirement (Rl), F ra d — F eq (T eq , V) for Tl = Tr = T eq . A related discussion is given in section Wl 

Note that, with only the relation JSJ, an arbitrary function of Tl and Tr can be added to the definition of F ra( i 
as F ra d = — (1/3) trP V + K(Tl,Tr). However F ra d should be an extensive variable due to the requirement (R3). 
Consequently, the arbitrary function disappears, K = 0. 



H. Chemical potential 



The chemical potential in general can be interpreted as a work needed to add a particle to a system under consid- 
eration. Because a photon is a collisionless particle, no work is required to add a new photon into the radiation field 
in the cylinder. This is the case either for equilibrium states or for steady states. Indeed, the chemical potential of a 
radiation field in equilibrium is zero. Therefore it is appropriate to set the SST chemical potential of a radiation field 
also as zero. 



I. Intensive nonequilibrium order parameter 

As a candidate for a nonequilibrium order parameter of a steady state, we can recommend the energy flow J, its 
areal density j = a (Tf; — Tr) and the temperature difference, 

r = T L - T R , (8) 

where r is defined as a positive quantity due to the assumption Tl > Tr . 

Note that, as mentioned at subsection IIII Al J is not a state variable. Therefore J is excluded from a state variable 
describing a nonequilibrium order. Next look at the energy flux j. If j is a consistent state variable, its conjugate 
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variable should also be denned in a consistent way. However the following calculation shows that the conjugate 
variable to j disappears for any value of j, 

dF rad _ 8a f 3 dT L dT R \ 8a f 1 1 \ 



This means that the conjugate variable to j can not be defined in a consistent way, and that j should also be excluded 
from a state variable describing a noncquilibrium order. 

On the other hand, r is obviously intensive and satisfies the requirement (R2), r — for Tl — Tr. Thus we adopt 
r as the intensive state variable for a nonequilibrium order of steady states. 

J. Temperature and extensive nonequilibrium order parameter 

We proceed to search for a definition of an SST temperature, T ra d- One may naively try to use the relation, 
dS ra d/ dE ra( i — 1/T ra d, to obtain the form of T ra d- In ordinary equilibrium thermodynamics for a radiation field, 
this partial derivative is calculated by fixing the extensive variable V. For our system, an extensive state variable of 
a nonequilibrium order of steady states should also be one of the fixed variables in the partial derivative. However, 
because such an extensive variable has not been specified yet, we can not evaluate the partial derivative here. Therefore 
we try to search for the form of T ra d with the other relation, 



9F ra d 



^rad ■ 

V,T 



v dT rad , 

Further we require that the temperatures Tl and Tr depend on T ra d and r, 

Tl = Tl (T ra d, t) , Tr = T R (T ra d, t ) . 
Then using the form of S ra d and F ra d we find, 

dF rad _ _ dT L dT R 

Therefore, with the relation r = Tl — Tr, we can set 

Tl = T r ad + a r + L(r) , Tr = T rad - (1 - a) r + L(t) , 

where a is a constant and L(t) is a function of r determined below. These temperatures give the form of T ra d as 

T rad = (l-a)T L +aT R -L(T). (10) 

For later use to determine the form of L(t), wc define an extensive noncquilibrium order parameter, as 

^_9F™, = 8a 3 3) 

or 6 c ' 

where L'(t) = di(r)/dr. The requirement (R2) gives the relation, 

V(T L = T R ) = £'(0) = i-a. 

Further the Legendre transformation should hold due to the requirement (R5), 

F r ad = F ra d{Trad, T, V) + T ra d S r ad + T * , (11) 

where F ra d is treated as a function of T ra d, t and V. Here note that, since r has the dimension of temperature and 
^ has that of entropy, the signature of the term r ^ in equation 111( 1 is positive. Substituting ^ and T ra d of (|ll)|) into 
(|ll[l . wc obtain 

L'T L -L'T R = L => L'(t) = ^ => L(t) = \t, 
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where A is a constant. Consequently we find a + A = 1/2, and T ra d of the form 110(1 becomes 

T rad = ~ {T L +T R ) . (12) 
The extensive nonequilibrium order parameter becomes 

* = ^ ( r £ 8 - t£ ) y = 1 ( s eg (T L , vo - s e9 (T fl , v) ) . (13) 

The form of T ra( 2 denotes that the third law holds in the SST for a radiation field if the third law of the equilibrium 
thermodynamics is retained for the left and right lids. 

K. Consistency of thermodynamic functions and 1st law 

From above results, it can be checked easily that the requirement (R4) is satisfied, 

8 2 F ^ B 2 F 8 2 F a 

-gy2~ = 0(> 0) : convex / -^ T - < , < : concave. 

Further, we find that the requirement (R6) is satisfied, 

~S rad dT rad + V dP rad - * rfr = , (14) 
Erad = T r ad S r ad - PradV + T * , (15) 

where P ra d is given by the SST pressure tensor P as equation ©. These are the Gibbs-Duhem relations. Then we 
obtain the first law, 

dErad = Trad dS r ad ~ Prad dV + T dUf . (16) 

Here note that, for example for an ordinary continuum system, the work term in the first law is given like Pij detj, 
where Pij is the stress tensor (which corresponds to the pressure tensor) and e,j is the strain tensor. The tensor e,j 
can be decomposed, = (1/3) tr(e) + e^J , where the trace tr(e) is the bulk strain and e^J is the constant volume 
strain, tr(e( c )) = 0. That is, for an ordinary continuum system, the work term becomes 

Pij deij = ^trPdV + Pij de[f , 

where dV = d[tr(e)] means the change of volume. With this observation, we can discuss the first law (|16|) as follows. 
In the context of the cylinder case, the constant volume strain means the change of the cross section A and the length 
L of the cylinder as A — > A A and L — > L/X, which denotes V does not change. Then, since all of the state variables 
obtained in this section HTll are characterised by Tj,, Tr and V, the constant volume strain can not cause a change 
of state variables like E ra d- That is, no interaction among photons in a radiation field causes zero energy change in 
the radiation field by any deformation of the system so long as the whole volume of the radiation field V does not 
change. Hence it is reasonable that the constant volume strain does not contribute to the first law l|ltj|) and only the 
bulk strain appears in the work term as —P ra ddV in (|16fl . The absence of the work term of the constant volume 
strain should be understood as a unique property of a radiation field. More discussions are given in section Ivl 
Finally, it follows that 



dSrad \ 1 
'V* T ra 



dE. 



ra 



This denotes that our definition of T ra d, S ra d, E ra d an d \£ are theoretically reasonable. 

As seen so far, the thermodynamic formulation (0th, 1st, 2nd and 3rd laws) can be obtained for the steady states 
of a radiation field in the cylinder case. 
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FIG. 5: The region B is excluded from the volume of a radiation field in local steady states. 



IV. CAVITY CASE 



A. Local steady state 



The cavity case (right one in figure ^) is treated in this section. The geometrical shapes of the inner black body 
and the cavity are arbitrarily set. For example, consider a case that the inner black body is a thin board, and count 
the number of photons coming from the board to a spatial point. Such a photon number at a point in front of the 
surface of the board is larger than that at a point on the plane of the board. Thus, in general, it is recognised that 
the number of photons coming from the inner black body to a point varies from point to point. The same is true 
of the number of photons coming from the outer black body. This means that the distribution function of photons 
is anisotropic and inhomogeneous as shown in next subsection. On the other hand, for example, the anisotropy and 
inhomogcncity of the velocity field of an ordinary fluid system have been traditionally handled well with the idea of 
local equilibrium in the ordinary fluid dynamics. Hence, after such a success of the idea of locality for the anisotropic 
and inhomogeneous system, we adopt the idea of local steady states. The radiation field in a sufficiently small region 
is in a steady state, but the steady state in one small region may be different from that in the other small region. 
Therefore the state variable in the cavity case should be defined as a function of a spatial point x, where the extensive 
variable is to be understood as a density. 

The idea of local steady states is not unique for a radiation field, but it also appears in formulating the steady state 
thermodynamics for dissipative systems!]]]. It is the inhomogencity of the temperature of the matter that raises the 
idea of local steady states in reference [tJ- On the other hand, it has already been revealed by reference [!j that, for 
a steady state of a hard core molecule gas system, the inhomogeneity of the temperature causes the inhomogeneity 
of the distribution function. Therefore the origin of the idea of local steady states in reference J7J seems to be the 
same as that in our radiation system. Further in the extended irreversible thermodynamics [3j, the idea of local 
nonequilibrium state is suggested, since, for example, some state variables are treated as functions of spatial point. 
Therefore, the idea of locality seems to be universal without respect to the nature of interactions among components 
of the nonequilibrium system under consideration, as though the formalism presented in this paper, for example the 
exclusion of the energy flux from state variables, is unique to steady states of a radiation field. 

Before proceeding to the search for definitions of state variables, let us note about the spatial volume filled with a 
radiation field. Because photons are collisionless, the region B shown in figure [S] should be excluded from the volume, 
V ra d, of a radiation field which is in local steady states. The radiation field in the region B is in an equilibrium state 
of temperature T out . 



Let pi n and p out denote the momentums of photons emitted at the inner black body and at the outer black body 
respectively. The directions from which the photons of pi„ and p ut can come to a spatial point x vary from point to 
point. Then the distribution function in the cylinder case Q is extended to the cavity case as 



B. Distribution function 



d{x,p) 



1 



(17) 



exp[huj/T(x, p)] — 1 ' 



where uj — pc/h, and T(x,p) is given by 




out 



for p = pi n at x 
for p = p out at x 
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FIG. 6: The geometrical factor, gi„(x), is shown. Another one is g ut(x) = 1 — gi„{x). 



C. Internal energy density 

The SST internal energy density, e ra d(x), is defined as the energy of photons at a point x, 



e r ad(x) = 2 J j^- huj d(x,p) = ^ J dQ xp T ln + J dQ xp T c 



t4 



where the factor 2 in the first line is due to the helicity of a photon, d£l xp is the solid-angle element in p-space at a 
spatial point x. Because the solid-angle of the inner black body measured from a spatial point x takes on different 
values according to the position of x, the numerical values of the integrals in the last line vary from point to point. 
Therefore we introduce the geometrical factors, gi n (x) and g ou t (x) , as 



9in{x) 



1 

47T 



<ifL 



, 9out(x) 



1 

47T 



(18) 



One of them, gi n (x), is shown in figure H3 These geometrical factors satisfy a relation, gi n {x) + g ut(x) = 1, and are 
positive, gin(x) > and g ut(x) > 0. Then the SST internal energy density becomes 



4cr 



e ra d(x) = — (gin(x) T^ n + g out (x) T* ut ) = g in (x) e eq (T m ) + g out (x) e eq (T out ) , 



(19) 



where e eq {T) = (4<j/c)T 4 is the equilibrium internal energy density of temperature T. This e ra d{x) satisfies the 
requirement (Rl), e rac i — e eq (T eq ) for Tj„ = T out = T eq . Further due to the requirement (R7), a total internal energy 
of a radiation field, E ra d, is given by a volume integral of e ra d{x), 

dx G-radix*) , 

d 

where the region like B in figure [S] should be excluded from the integral region. It is obvious that the cylinder case is 
recovered by setting g in {x) = g ut{x) = 1/2. 



D. Entropy density 

As for the cylinder case, we try referring |14| |Tfij and §55 in to define the SST entropy density, s ra d(x), as 
s r ad(x) = 2 J -p- [(1 + d(x,p) ) ln(l + d(£,p) ) — d(x,p) lnd(x,p)] 

= (9in(x) T t l + g ut(x) T^ ut ) = g in (x) s eq (T in ) + g out (x) s eq (T out ) , (20) 

where s eq (T) — (16(t/3c) T 3 is the equilibrium entropy density. This s ra d{x) satisfies the requirement (Rl), s ra d(x) = 
s e q{T eq ) for Ti n = T out = T eq . Due to the requirement (R7), a total entropy of a radiation field, S ra d, is given by a 
volume integral of s ra d{x), 



Srad — / dx S ra d(x) . 
JVrad 

The cylinder case is recovered by setting gi n (x) = g ou t{x) = 1/2. As discussed at the last paragraph in subsection 
IIII Dl s ra d(x) is just one candidate for a well defined SST entropy density. Hence, to be exact, we observe whether 
we can construct a candidate for a consistent SST for a radiation field with s ra d{x) of equation l|2U|l . 
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E. 2nd law with a relaxation process 

This subsection shows that s ra d{x) of equation l|20l) satisfies the requirement (R8). We set the outer black body 
being covered with a heat insulator and consider that each black body passes a sequence of equilibrium states and 
the radiation field passes a sequence of steady states during the relaxation process of the whole system. Due to the 
conditions (CI) and (C3), the ordinary thermodynamic relations hold for the black bodies, and the volumes of the 
two black bodies does not change. Further due to the requirement (R7), the total energy E ra d and the total entropy 
S r ad are given by E tot = E in + E out + E rad and S toi = S in + S out + S rad respectively, where E in and S in are the 
equilibrium internal energy and entropy of the inner black body, and E out and S out are those of the outer black body. 

Here it should be noted that the inner black body may move due to the radiation pressure. However, because of 
the condition (C3), the radiation volume V ra d does not change, dV ra d = 0. Further it is appropriate to assume that 
the velocity of the inner black body is very small in comparison with the speed of light. Hence the special relativistic 
Doppler effect and the kinetic energy of the inner black body can be ignored in the analysis of the relaxation process. 
And it should also be pointed out that the inner black body rests somewhere in the cavity when the whole system is 
in a total equilibrium state. 

With the aid of the energy conservation and dV ra d — 0, a calculation similar to the one done in subsection IIII El 
gives the following inequality, 

dStot = tp — dTi n + — — dT ou t H ( Gi n T^ n dTi n + G ou t T 2 ut dT 0U t ) 

in 1 out C 

= ( tp Tf— J ( Cout H G out T 3 ut ) dT out > , 

V 1 out J-in J \ C / 

where Gi n — j v ^ dx 3 gin^x), G ou t = f v dx g ut{x), and Ci n and C ou t are the heat capacities of the inner and outer 
black bodies respectively. The inequality, dStot > 0, holds for both cases Ti n > T ou t and Ti„ < T ou t, and the equality, 
dStot = 0, holds for the equilibrium case Ti„ = T ouf . This indicates that s ra d{x) of equation (|20|l is a candidate for a 
well defined SST entropy. 



F. Pressure 



One may naively expect that the pressure of a steady state is a global quantity, since the pressure gradient in an 
ordinary dissipative system accelerates the components of the system to cause a dynamical evolution in the system. 
However, since no dissipation exists among photons, this discussion can not be applied to a radiation field. As seen 
below, the spatial dependence of d(x, p) makes the SST pressure tensor for the cavity case a function of x. 

As for the cylinder case, we define the SST pressure by the tensor, 

Pn„(x)=2j^T- c CuC„d(x,p). (21) 

It can be checked by the same calculation in subsection IIII Fl that this P(x) satisfies the requirement (Rl). 



G. Free energy density 

The discussion of the free energy given in subsection IIII Gl is also applied for the cavity case. Therefore, since we 
are searching for a density of an extensive variable, the SST free energy density, f ra d{x), is defined as 



frad (x) = - i trP (x ) = - ^ J Huj d(x , p) = - i e rad 



(x) 



Hence we obtain 



frad{x) = --^ (gin(x) T£ + g ut{x) T D \ t ) = g in {x) feq(T in ) + g ut{x) f eq {T out ) , (22) 
o c 

where f eq (T) = -(4cr/3c)T 4 is the equilibr ium free energy density. This f ra d(x) satisfies the requirement (Rl), 
frad = feq(T eq ) for T in = T ou t = T eq . The cylinder case is recovered by setting g m {x) = g ou t(x) = 1/2 and integrating 

f radix) Over V rad - 
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H. Chemical potential 

The discussion of the chemical potential given at subsection IIII HI is also true of the cavity case. Hence we set the 
SST chemical potential to be zero in the cavity case, too. 



I. Intensive nonequilibrium order parameter 



There are three candidates for a nonequilibrium order parameter; the temperature difference r = T in — T out , the 
energy flow J — a (T in — T^ ut ) Ai n , where Ai n is the surface area of the inner black body, and the third candidate 
is the areal density of J at the surface of the inner black body ji n = a (T in — T^ ut ). If J and jin are adopted as 
state variables, it is impossible to recover the cylinder case by setting gi n (x) = g ou tix) = 1/2, since J and ji n in the 
cylinder case are not state variables. Therefore it is also appropriate for the cavity case to adopt r as a state variable 
of a nonequilibrium order. Here r is intensive in the cavity case, too. More discussion on the definition of r is given 
in section Ivl 



J. Temperature and extensive nonequilibrium order parameter 



To find a definition of the SST temperature, T rad (x), we follow the same discussion and procedure given in subsection 
IIII .11 The same calculations until (|10|l hold for the cavity case except for that the constant a and the function L(t) 
are modified to depend on a spatial point, a{x) and L(x; r). Then we define an extensive state variable, ip(x), which 
is conjugate to r, 

m = - = ^ ( ( a + L > ) gin T? + ( -(1 - a) + L' ) g out T a l t ) , 

where L' = 9L(x;t)/9t. The requirement (R2) gives the relation, 

t/j(x) = for T in = T ou t => L'(x; 0) = g ou t(x) - a(x) , 

where gi n + g out — 1 is used. Further due to the requirement (R5), the Legendre transformation (|11|) should also hold 
in the cavity case, 

e ra d{x) = frad(x; T rad (x),T) + T rad (x) s rad {x) +Tip(x) . (23) 
Substituting an d T rad of the form l|10|l into Ij23(l . we obtain 

L' T in — L' T out = L => L(x;t) = X(x)t , 
where A is an unknown function. Consequently we find a(x) + X(x) = g ut(x), and T rad of form (|10fl becomes 

g in {x) T in + g ut(x) T ut ■ (24) 
The nonequilibrium order parameter ip(^) becomes 

i>{x) = ^7 9in(£) 9out{%) ( T il- T out) = 9in{x) g out {x) (s eq (T m ) - s eq (T ou t)) ■ (25) 

The cylinder case is recovered by setting gi n (x) = g ut{x) — 1/2 and integrating ip(x) over V rad - The form of T rad (x) 
shows that the third law holds in the SST for a radiation field if the third law of the equilibrium thermodynamics is 
retained for the inner and outer black bodies. 



K. Consistency of thermodynamic functions and 1st law 

From the above results, the requirement (R4) is automatically satisfied, 

d 2 frad(x) d 2 f rad (x) 
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Further the requirement (R6) is satisfied, 



-s ra d(x) dT rad (x) + dP rad {x) - ip(x) dr | fcfixcd = , (26) 
e ra d{x) = T rad (x) s rad (x) - P ra d(x) + t ip(x) , (27) 



where 



Prad(x) = i trP(x) = ^- (g ln (x) + g ut{x) T^ ut ) . 
Then the first law is obtained, 

de rad (x) | fcfixod = T rad (x) ds rad (x) +Tdip(x) | 5=fixcd . (28) 

Finally it follows 

'ds r ad(x)\ 1 



de ra d{x) J ^fflg T rad (x) 

This indicates that T rad (x), s rad (x), e rad (x) and ip(x) is theoretically reasonable. 

As seen so far, the thermodynamic formulation (Oth, 1st, 2nd and 3rd laws) can be extended to the steady states 
of a radiation field in the cavity case. 



V. SUMMARY AND DISCUSSIONS 



Summary: Based on the collisionless nature of photons, we have been able to find the distribution function of 
photons in a two-temperature steady state. Consequently, using the Shannon- type entropy JSJ, a candidate for a 
consistent two-temperature steady state thermodynamics (SST) for a radiation field in vacuum has been constructed 
in the cylinder case, where state variables have been defined globally. It was impossible in the cavity case to define 
state variables globally, since the system is anisotropic and inhomogeneous. However, after introducing the idea of 
local steady states which is motivated by the anisotropy and inhomogcncity of the system, we have been able to 
extend the SST to the cavity case. As mentioned in the second paragraph in subsection IIII Bl the state variables in 
the SST for a radiation field are given by an appropriate linear combination of two equilibrium values. 

It should be emphasised that, contrary to the cases of ordinary dissipative systems 0> the energy flux can 
not work as a consistent state variable for a radiation field 0, as has been explicitly shown in the equation 10. 
Concerning this issue, we have found in this paper that the consistent state variables of the nonequilibrium order of 
steady states are the difference of the temperatures of two black bodies r and its conjugate variable ip(x). A somewhat 
detailed discussion of the energy flux is put in appendix ^ where the reason why an energy flux can not behave as a 
state variable is explored. 

Special properties of a radiation field: Other than the exclusion of the energy flux from consistent state 
variables, there are other special properties of a radiation field. As mentioned below the first law (|16f) . the constant 
volume strain has no physical meaning in the first law for a radiation field. Further the first law (|16fl together with 
the Legendre transformation l|llf> gives, 

dFrad = -S rad dT rad - P rad dV - * dr . 

This shows that the components of the SST pressure tensor P can not be obtained from the SST free energy F rad . 
Such a limitation on the information obtained from the free energy should be understood as a unique property of a 
radiation field. 

However, if one proceeds to a general relativistic extension, it may be possible to obtain the components of P from 
the free energy. The pressure tensor should be extended to the so-called stress- energy-momentum tensor T in the 
context of general relativity. The tensor T is related to the deformation R (curvature tensor of rank two) of the 
spacetime via the Einstein equation, and the tensor R can be interpreted as the strain of the spacetime raised by the 
radiation field. That is, the work term may be given as T : c£R. Hence we may obtain T by the partial derivative of 
the free energy by R. 

The distance of steady states from equilibrium states: Let us remark on the method of finding the equilibrium 
limit and on the definition of the intensive nonequilibrium order parameter t. In constructing the SST for the cavity 
case, we have considered an equilibrium limit as Ti„ — T out — > 0. However there is another method of finding the 
equilibrium limit which is to remove the inner black body. This limit is taken with gi n (x) — > (<=> g ut{x) — ► 1). With 
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this limit, r remains at a constant value Ti n — T out before gi n (x) reaches zero, but r becomes T out discontinuously at 
gin{x) = 0. That is, r does not disappear at this equilibrium limit, and not satisfy the requirement (R2). However, 
because i/j(x) ~~ * f° r both equilibrium limits Ti„ — T out — > and gi n (x) — > 0, the violation of the requirement (R2) 
by r is harmless and our formulation of the SST remains consistent. 

If one wishes to modify the definition of r to satisfy the limit r — > as gi n (x) — > 0, one method is to re-define 
as t = (gi n gout) q (Ti n — T ou t), where q is a constant and < q < 1. Consequently ^(x) should be modified to 
V'(x) = (gingout) 1-11 (s e q(Tin) — s eq (T out )). However as discussed in the previous paragraph, q — does not cause 
inconsistency in formulating the SST. Therefore we adopt q = as the simplest definition of t. 

The nonequilibrium order parameters r and i/)(x) seem to give some measure of distance of steady states from 
equilibrium states. It is appropriate to require that the distance of steady states from equilibrium states becomes 
zero with both equilibrium limits, T in — T out — > and gi n (x) — > 0. Note that, when we adopt the simplest definition 
of r = Ti n — T out , this t can not disappear with the limit gi n {x) — > 0. Therefore it is plausible to consider ip(x) 
as the measure of distance of steady states from equilibrium states, where this distance should be interpreted as a 
kind of averaged distance from two equilibrium states of temperatures Tj„ and T ou t- The larger the value of ip(x), 
the more distant a steady state from two equilibrium states. In fixing the values of T n and T ou t, the distance of 
a steady state from two equilibrium states can vary by changing the geometry of system. Here note an inequality, 
< a (1 — a) < 1/4 for < a < 1, where the maximum value 1/4 is taken with a — 1/2. Then we find that ip(x) 
takes the largest value for the case gi n (x) = g ut(x) = 1/2. On the other hand, the cylinder case is recovered from the 
cavity case by setting <^ n (x) = g ut{x) — 1/2. Consequently, it is concluded that the radiation field in the cylinder 
case is in the most distant steady state from two equilibrium states. This implies that the higher the geometrical 
symmetry of the system, the more distant the steady state from two equilibrium states. 

Future improvements: We have assumed the condition (CI). In a practical experiment, the condition (CI) 
may imply that the time scales of thermalisation of the two black bodies are very short. Then a heat flux probably 
penetrates well into the black bodies, and their surfaces are affected by the heat flux. In such a case, the distribution 
functions and (|17|l should be extended to nonequilibrium ones. Therefore, because the radiation field is completely 
characterised by the temperatures at the surfaces of the two bodies, the nonequilibrium effect of the surfaces probably 
becomes essential for the SST for a radiation field. 

If we consider T ou t — for the cavity case, then no energy flux comes into the inner black body and no nonequilibrium 
effect on its surface arises. Such a case can be considered in a general relativistic system. It has already been known 
that a black hole emits a thermal radiation, that is, a black hole can be considered as a black body 17] 18] . Then, if 
a black hole is isolated from other celestial bodies, no energy flux may come into the black hole. Therefore the cavity 
case with T ou t = can be applied to the system composed of a black hole and its thermal radiation. However before 
proceeding to such an application, we may need to carry out a relativistic extension on the SST for a radiation field. 

Applications: In laboratory experiments which need to keep an instrument isolated, a method used frequently for 
heat insulation is to keep the instrument be surrounded by vacuum, where the instrument corresponds to the inner 
black body in the cavity case. A total value of the nonequilibrium order of a radiation field in the surrounding vacuum 
is estimated as ^ = f v dx 3 ip(x) ~ (T in — T Q \ t ) V ra d x 1CT 8 erg/K, where we set g in ~ g ut ~ 0(1). For example, 
consider the case that the instrument is sufficiently near an equilibrium state whose equation of state is of an ideal 
gas type. Then the entropy of the instrument is estimated as Si n ~ (E in + PinVin)/T in ~ uNa ks ~ n x 10 7 erg/K, 
where n is a numerical factor which denotes the size of the instrument and Na is the Avogadro's number. Then we 
find, *S? / 'Sin ~ (T in ~ T Q 3 ut ) V r ad/n x 1CT 15 . This denotes that an effect of a radiation field in the surrounding vacuum 
does not come into account unless we have sufficiently large difference of temperatures, r > 10 5 K, or sufficiently large 
volume of the surrounding vacuum, V ra d ^ 10 15 cm 3 , or sufficiently small instrument, n < 10~ 15 . That is, the heat 
insulation with surrounding an instrument by vacuum is very effective in laboratory experiments. 

A situation in which an effect of the SST for a radiation field becomes important can be considered in astrophysics 
and cosmology, since the volume of vacuum space is sufficiently large in inter-stellar spaces and temperatures of 
celestial bodies can be sufficiently large. When a radiation field in an inter-stellar space is considered, the condition 
(C2) is violated and we need to take the effects of the retarded time on photons into account. However, when a 
radiation field inside a star is considered, it is possible to obtain a notable implication for astrophysics with the 
present form of the SST for a radiation field. A somewhat detail explanation of this astrophysical implication is 
summarised in appendix [BJ The point of this implication is that, in comparison with the usual analysis of a self- 
gravitating gas system using the ordinary (local) equilibrium thermodynamics, the SST for a radiation field tends 
to prevent the so-called gravothermal catastrophe from occurring. The gravothermal catastrophe is one of the most 
significant phenomena of self-gravitating systems, which is the essential physical basis to characterise the formation 
and evolution of astrophysical objects like stars and galaxies [l| |l6j. Therefore the SST for a radiation field may give 
some useful corrections or developments in the theories of stellar structure and other astrophysical systems. 
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APPENDIX A: ENERGY FLUX 



Let us consider the cavity case. Look at a spatial point x in the radiation field, and take the sum p to t of the 
momentums pi n and p out of all photons passing the point x at one moment. The direction of a net energy flow 
through x at that moment is parallel to ptot- Then we can define the energy flux vector as 

3 =jft, 

where ft is the unit vector parallel to ptot at x, and the magnitude of this vector j is the energy flux given by 

j = 2j -j^- h lu d(x , p) c cos 4> , 

where the factor 2 is due to the helicity of a photon and <fi is the angle between p and ft given by p cos <fi = ft- p. With 
the relation / °° dxx 3 /(e x — 1) = n 4 /15, we obtain 

j(x) = cr ( 7m (f) T£ + j out (x) T* ut ) , 

where "fm{x) and 7 ott t(x) are the geometrical factors given as 

Jm (x) = - [ dtt xp cos (f) , 7ottt(a?) = — / dtt xp cos 4>, 

where dfl xp is the solid-angle element in p-space at x. These satisfy the relation, 7,„(f) + Jout(x) = 0, then j = at 
the equilibrium limit T in — T out . The factors 7i„(af) and "f ut{%) are obviously different from the factors gi n (x) and 
g ou t{x) which appear in the form of the state variables given in section ITVl This difference may be interpreted as one 
reason why the energy flux does not work as a consistent state variable. 



APPENDIX B: IMPLICATION FOR ASTROPHYSICS 



A remarkable property of a self-gravitating system is the negative heat capacity of a region around the centre of 
mass due to the virial theorem and the absence of the lower bound of the gravitational potential (oc — 1/r, where r 
is the distance from the centre of mass) pj . The significant phenomenon caused by the negative heat capacity is 
the gravothermal catastrophe, which is, for example, the essential basis to cause and maintain the nuclear reactions 
in the core of the star. The simplest toy model of the gravothermal catastrophe is given by setting the volume of 
the radiation field in our cavity case as zero (see figure EJ) , and setting the heat capacity of the inner black body 
as negative, Ci n < 0. Let us consider the case T^ n > T out . Then, while the temperature of the outer black body 
continues to increase due to the incoming heat from the inner black body, the temperature of the inner black body 
does also continue to increase. Further if the relation \Ci n \ < C ou t holds, the rate of temperature increase of the inner 
black body is faster than that of the outer black body. Then the difference of the temperatures continues to increase 
without upper limit. This is a rough illustration of the gravothermal catastrophe. For an example of the case of a star 
formation, the inner black body corresponds to the core of the star and the outer black body to the other part of the 
star, and the temperature increase of the core due to the gravothermal catastrophe triggers the nuclear interactions 
in the core. Here note that, since the virial theorem is used for deriving the negative heat capacity, the gravothermal 
catastrophe is a consequence of a local equilibrium treatment of the system under consideration. And, even if the 
radiative energy transfer is also included in the above discussion, the gravothermal catastrophe should be concluded 
whenever the radiative energy transfer is treated in a traditional manner mentioned as in the beginning of section [I] 

Let us again place a radiation field between the inner and outer black bodies, and examine the relaxation process 
with the negative heat capacity of the inner black body Cm < 0. The same calculation given in subsection II V El holds 



18 




FIG. 7: A toy model for the gravothermal catastrophe. The heat capacity of the inner black body is negative. 



in this case. Here look at the total energy and its differential, 



— | E-,„ H — — Gi n T^ \ + ( E ou t H — —G utT^ ut j . 
dEtot = ( Ci n H Gi n T.?. I dTi n + I C ou t H G ou t T^.f I dT D 



tot — '-'in i '-'in J- in I UJ in T I <-^cmt i "-^out -'-out I 111 out i 

V c / V c / 

where the terms including Gi„ or G ou t come from the SST internal energy E ra d- This enables us to consider our 
system as being composed of two subsystems X and Y whose heat capacities Cx and Cy are given as 

Cx = Ci n H Gi n T in , Cy — C ou t H G ou t T o t , 
c c 

and the temperatures of the subsystems are Ti n for X and T out for Y. For simplicity set Cj n = constant < and 
> T out . In this case, energy flows from the subsystem X to Y, and the heat capacity Cx of X is negative for 
a sufficiently low temperature of the inner black body T in . Therefore T in increases. However, after T in increases 
sufficiently, the heat capacity Cx turns from negative to positive. This denotes that the evolution of Ti n during the 
relaxation stabilises at a critical temperature T* n determined by Cx = 0. Hence it is expected that the whole system 
settles in the equilibrium state of the critical temperature T* n . This discussion is applicable to the case |Cj„| < C out , 
which results in the gravothermal catastrophe when the system under consideration is treated with a local equilibrium 
theory. Therefore it is implied that the SST for a radiation field tends to prevent the gravothermal catastrophe from 
occurring. The theory of the evolution of the stellar structure and the rate of star-formation in a galaxy may be 
corrected if the SST for a radiation field will be applied to the radiation field inside a star. 
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